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Abstract — In cognitive radio (CR) networks with multiple- 
input multiple-output (MIMO) links, secondary users (SUs) 
can exploit "spectrum holes" in the space domain to access 
the spectrum allocated to a primary system. However, they 
need to suppress the interference caused to primary users 
(PUs), as the secondary system should be transparent to the 
primary system. In this paper, we study the optimal secondary- 
link beamforming pattern that balances between the SU's 
throughput and the interference it causes to PUs. In particular, 
we aim to maximize the throughput of the SU, while keeping 
the interference temperature at the primary receivers below a 
certain threshold. 

Unlike traditional MIMO systems, SUs may not have 
the luxury of knowing the channel state information (CSI) 
on the links to PUs. This presents a key challenge for a 
secondary transmitter to steer interference away from pri- 
mary receivers. In this paper, we consider three scenarios, 
namely when the secondary transmitter has complete, par- 
tial, or no knowledge about the channels to the primary 
receivers. In particular, when complete CSI is not available, 
the interference-temperature constraints are to be satisfied 
with high probability, thus resulting in chance constraints that 
are typically hard to deal with. Our contribution is fourfold. 
First, by analyzing the distributional characteristics of MIMO 
channels, we propose a unified homogeneous quadratically 
constrained quadratic program (QCQP) formulation that can 
be applied to all three scenarios, in which different levels of 
CSI knowledge give rise to either deterministic or probabilistic 
interference-temperature constraints. The homogeneous QCQP 
formulation, though non-convex, is amenable to semidefinite 
programming (SDP) relaxation methods. Secondly, we show 
that the SDP relaxation admits no gap when the number 
of primary links is no larger than two. A polynomial-time 
algorithm is presented to compute the optimal solution to the 
QCQP problem efficiently. Thirdly, we propose a randomized 
polynomial-time algorithm for constructing a near-optimal 
solution to the QCQP problem when there are more than 
two primary links. Finally, we show that when the secondary 
transmitter has no CSI on the links to primary receivers, the 
optimal solution to the QCQP problem can be found by a 
simple matrix eigenvalue-eigenvector computation, which can 
be done much more efficiently than solving the QCQP directly. 
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I. Introduction 

Cognitive radio (CR), which allows secondary users (SUs) 
to opportunistically utilize the frequency spectrum originally 
assigned to licensed primary users (PUs), is a promising 
approach to alleviate spectrum scarcity IH. In CR networks 
with single-antenna nodes, SUs can transmit only when it 
detects a spectrum hole in either time or frequency domain, 
so as to avoid causing harmful interference to PUs 121, O. 
Such schemes, however, only work when the primary system 
severely underutilizes the assigned spectrum. Otherwise, 
the secondary system would not have adequate chances to 
access the wireless channel. 

Recent development in multiple-input multiple-output 
(MIMO) antenna techniques opens up a new dimension, 
namely space, for co-channel users to coexist without 
causing severe interference to each other H. Indeed, in 
CR networks where stations are equipped with multiple 
antennae, SUs can transmit at the same time as the PUs 
through space-domain signal processing. The nature of CR 
networks gives rise to several challenging issues that do not 
exist in traditional MIMO systems. First, SUs are solely 
responsible for suppressing the interference they cause to 
PU receivers, as the primary system should not be aware of 
the existence of the secondary system. That is, we cannot 
rely on the PUs to do receiver-side interference cancellation. 
Secondly, SUs may not have the luxury of knowing the 
channel state information (CSI) on the links to PUs, as the 
primary system would not deliberately provide their channel 
estimation to the secondary system. This imposes difficulty 
on transmitter-side pre-interference cancellation at SU trans- 
mitters. It is therefore necessary to revisit space-domain 
signal processing in the context of MIMO CR networks. In 
particular, SUs need to configure their beamforming patterns 
in a way that balances between their own throughput and 
the interference they cause to PUs. 

Multi-antenna CR networks were recently studied in IS)- 
Q- Assuming that CSI on all links is perfectly known 
to the SUs, ID formulates the SU beamforming problem 
as a non-convex optimization problem. A semi-distributed 
algorithm is proposed to obtain a local optimal solution 
to the problem. On the other hand, Q assumes that the 
PU can act as a scheduler for SU transmissions. Under 
this idealistic assumption, an opportunistic orthogonaliza- 
tion scheme is proposed. In |[6j, Zhang and Liang studied 
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capacity-achieving transmit spatial spectrum for a single 
SU, assuming that the SU has full CSI and there is no 
interference from PUs to the SU. Insightful solution methods 
are proposed to provide better intuition that may not be 
obtainable from a numerical optimization perspective. The 
issue of imperfect CSI estimation is considered in ||8l- 
lITOl for multiple-input single-output (MISO) CR networks. 
Therein, robust optimization problems are formulated to 
ensure the service qualities for both SUs and PUs are 
satisfactory in the worst case. 

In this paper, we study the problem of optimal secondary- 
link beamforming. Specifically, we aim to maximize the 
throughput of the SU under the constraint that the inter- 
ference to PU receivers is below a certain threshold. In 
contrast to the previous work fSj-E], we consider three 
practical scenarios: (1) when the SU transmitter knows 
both the channel matrices from it to PU receivers and the 
beamforming patterns at PU receivers; (2) when the SU 
transmitter does not know the beamformer at PU receivers; 
(3) when the secondary transmitter knows neither the chan- 
nel matrices to PU receivers nor the beamforming patterns 
at PU receivers. Note that the deterministic interference- 
temperature constraints could be too stringent if the SU 
does not have full CSI (which is the case in the second and 
third scenarios), as the SU will need to consider the worst- 
case channel reaUzation when configuring the beamformer. 
Fortunately, many wireless applications can tolerate occa- 
sional dips in the service quality. In order to have a more 
efficient utilization of the spectrum, we can take advantage 
of this opportunity and replace the deterministic constraints 
by probabilistic interference constraints (also referred to 
as chance constraints). Chance constraints, however, are 
typically tougher to deal with than deterministic constraints. 
The contribution of this paper is fourfold: 

• We propose a unified homogeneous quadratically con- 
strained quadratic program (QCQP) formulation that 
can be applied to all three scenarios mentioned above. 
In particular, the homogeneous QCQP formulation 
can accommodate both deterministic and probabilistic 
interference-temperature constraints. 

• The QCQP formulation, though non-convex, is 
amenable to semidefinite programming (SDP) relax- 
ation methods. We show that the SDP relaxation admits 
no gap with the true optimal solution when the number 
of PUs is no larger than two. A polynomial-time 
algorithm is presented to compute the optimal solution 
efficiently. 

• When there are more than two PUs, we propose a ran- 
domized polynomial-time algorithm that can produce a 
provably near-optimal solution. Numerical results show 
that the solution produced by our algorithm almost 
achieves the optimal value. 

• In the third scenario where the SU transmitter knows 
neither the channel matrices nor the beamformer at 
PU receivers, we show that the optimal solution can 
be obtained very efficiently through a simple matrix 



eigenvalue-eigenvector computation. That is, there is 
no need to solve the QCQP problem in this case. 

We should emphasize that the incomplete CSI here is not 
to be confused with that in ISl- lfTOl . where it is assumed that 
the SU knows the CSI on all links, except that there may be 
uncertainty in the channel estimation. Robust optimization 
techniques are employed in these papers to deal with the 
worst-case channel realization. The case where only channel 
statistics is known to the SU transmitter is also studied in 
ifTOl . This is similar to a special case of the third scenario 
considered in our paper, when the receivers have only one 
antenna. 

The rest of this paper is organized as follows. The system 
model is described in Section|IIl In Section|IIIl we formulate 
the SU beamforming problem as a series of homogeneous 
QCQP problems for different scenarios. The SDP relaxations 
of these homogeneous QCQP problems are then introduced 
in Section HVl In Section[V] a polynomial-time algorithm for 
finding an optimal solution to the QCQP problem when the 
number of primary users is no larger than two is presented. 
In Section IVII we propose a randomized polynomial-time 
algorithm to find a near-optimal beamforming solution when 
there are more than two PUs. In Section IVIII the perfor- 
mance of the proposed schemes is evaluated via simulations. 
Finally, the paper is concluded in Section IVIIII 

II. System Model 

A. System Setup 

In this paper, we consider a CR network in which a 
secondary user intends to share the spectrum with a primary 
system consisting of K primary links. We shall discuss 
the possibility of extending the proposed approach to a 
multiple-secondary-link scenario in Section IVIIII In the 
sequel, we use the subscript 5* to denote the secondary 
link and the subscript k to denote the fc*'' primary link. 
Let Ms (or Ns) and Mk (or N^) denote the number of 
transmit (or receive) antennae of the secondary and primary 
links, respectively. We use Hs,s G ^JVsxMg j-q denote 
the channel matrix from the secondary transmitter to the 
secondai-y receiver, and H^^s e c^'^''^^^', Hs,fc € C^^''^'^ 
and H/jj e c^^x^^j to denote the channel matrices 
from the secondary transmitter to the fc*'' primary receiver, 
from the A:*'' primary transmitter to the secondary receiver, 
and from the j*'' primary transmitter to the fc*'' primary 
receiver, respectively. We assume a Rayleigh fading and rich 
scattering environment, so that the entries of the channel 
matrices are independently and identically distributed (i.i.d.) 
complex Gaussian random variables with zero mean and unit 
variance. As pointed out in ifTTIl . lfT2l . in an interference- 
limited environment, each active link should only transmit 
one data stream at a time to avoid excessive interference 
to other links. In this case, we use scalars xs and Xk to 
denote the transmitted signals by the secondary transmitter 
and the k*^ primary transmitter, respectively. Without loss 
of generality, suppose that E[|a;sp] ~ 1 and E[|a;fcp] = 1. 
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Let ts and rs denote the beamforming vectors at the 
secondary transmitter and receiver, respectively. Likewise, 
let tfc and be the beamforming vectors at the fc*'* primary 
transmitter and receiver In particular, we have ||t5||2 = Ps 
and ||tfc||| = Pk, where Ps and Pk are the transmit power 
of the secondary and the fc*'' primary links, respectively. 
Without loss of generality, we normalize the receive beam- 
forming vectors so that UrslH = 1 and ||rfc||| ~ 1. Then, 
the received signal at the secondary receiver after receive 
beamforming is 

K 

ys = ^/as^rfHs,stsa:s+X! V^sjrf Hs.fctfeXfc+rfns, 

fe=i 

where as.s and as.k denote the path losses from the 
secondary transmitter to the secondary receiver and from the 
k*^ primary transmitter to the secondary receiver, respec- 
tively, and ns ~ CJV{0, NqI) denotes a circular complex 
additive Gaussian noise vector at the secondary receiver As 
such, the signal to interference and noise ratio (SINR) on 
the secondary link is given by 



is = 
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The secondary link's transmission causes an interfer- 
ence signal ^afc.sr^Hfe^gtsxs at the output of the fc*'* 
primary receiver, resulting in an interference power of 
afc.sjr^Hfe^stsI , where a^.s is the path loss from the 
secondary transmitter to the fc*'' primary receiver 

Before leaving this subsection, we emphasize that channel 
matrices on different links are independent. Furthermore, the 
beamforming vectors and are solely determined by 
the channels between the nodes in the primary system, as 
the secondary system is transparent to the primary system. 
Therefore, and are independent of 5, H5 5, and 

B. Objective and Assumptions 

In this paper, we aim to find for the secondary link the 
optimal beamforming vectors ts and rs so that the SINR on 
the secondary Unk is maximized, while the interference to 
primary link k is below a tolerable threshold e^. Mathemat- 
ically, this can be formulated as the following optimization 
problem: 

max 75 (2a) 
S.t. afc.s|rfHfc,sts|' < efc = 1, . . . , Apb) 

||ts||2 < Ps,max, (2c) 

where Ps.max is the maximum transmission power of the 
secondary link. Throughout this paper, we assume that 
the path losses as.s^ cts.k, ctk.s, and akj do not vary 
significantly within the time period of interest, and hence 



are known to all stations. It is also reasonable to assume 
that the secondary user knows its own channel H5 5 at 
both transmitter and receiver sides. Moreover, the secondary 
receiver can estimate Hg fetfc for all k by overhearing the 
transmission of primary transmitters. 

In practice, however, the secondary transmitter may not 
know the CSI on the links between primary receivers, 
because primary receivers would not purposely provide CSI 
to the secondary system. In this paper, we are interested in 
the following three different scenarios: 
Scenario 1: The secondary transmitter has perfect knowl- 
edge of the vector gVk- 

This scenario corresponds to a time division duplex 
(TDD) system in which channels are reciprocal and the 
primary receivers use the same beamforming vectors for 
both reception and transmission. In this case, the secondary 
transmitter can estimate H^^rfc by overhearing the trans- 
mission of primary receivers. 

Scenario 2: The secondary transmitter knows H^- 5 but not 
rfc- 

This scenario corresponds to a TDD system in which 
primary receivers do not use the beamforming vector rfc 
for transmission. 

Scenario 3: The secondary transmitter has no knowledge 
about Hk.s and r^. 

This scenario corresponds to the case where the secondary 
link has no way to estimate the channel from the primary 
receiver. 

Note that in both Scenarios 2 and 3, constraints (l2b] i are 
no longer well defined. Indeed, for any given ts, the value of 
the left-hand-side is uncertain to the secondary transmitter, 
as the realizations of and/or Hk.s are unknown. There- 
fore, a revision is necessary. One way is to guarantee that the 
constraints are always satisfied regardless of the realizations 
of Vk and Hfe 5. Then, ( l2b] i can be replaced by 

max Q;fc,s I rf Hfe^sts I ^ < efe \/k^l,...,K, (3) 

where the maximization is taken over for Scenario 2, and 
over Tfe and 'H.k,s for Scenario 3. 

Besides the worst-case guarantee, in practical applica- 
tions, we may allow the interference to exceed a certain 
threshold Ck with a small outage probability Sk- In this case, 
(l2bl ) can be replaced by 

Pr{afe,s|rfHfe,sts|' < e^} > 1-4 yk^l,...,K, 

(4) 

where the probability is taken over r^; for Scenario 2, and 
over rk and Hfc 5 for Scenario 3. Note that (|4|i is equivalent 
to (O when Sk — for all k. 

C. Distribution of 

In this subsection, we discuss the probability distribution 
of the beamforming vector at PU receivers. The results 
will be useful later when we address the chance constraints 
in (Hli. 
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As mentioned, is solely determined by the channels 
between primary nodes. Consider the matrix 

H = [-y/afc,iHfe^iti, . . . , ^afc^fcHfc^fctfe, 

:= [hi,...,hfe,...,hif], 

where akj is the path loss from the j*'* PU transmitter to 
the fc*'' PU receiver, and hi,...,hx are the columns of 
H, which are independent of each other Let H(_fc) be the 
Nk X {K — 1) matrix obtained by deleting the k*^ column 
of H. Then, in general, the vector takes the form 

rk^f3kWkhk, (5) 
where Bi, = -- — ,, is a normalization factor that ensures 

^'^ ||Wfchfc||2 

Ikfclli — 1' ™d Wfe is a random Hermitian matrix that is 
independent of hfe. In particular, we have = I for the 
matched-filter (MF) receiver, 

W, = I - H(_,)(Hf_,)H(.,))-iHf_^ (6) 

for the zero-forcing (ZF) receiver, and 

Wfc = (H(_,.)Hf_fe)+iVoI)-i (7) 

for the minimum-mean-squared-error (MMSE) receiver. Be- 
fore proceeding, let us state a definition and introduce two 
assumptions. 

Definition 1. A random vector x is called a normalized 
complex Gaussian vector if :x. = where z ^ CM{Q, I). 

In particular, ||x||2 = 1. A normalized complex Gaussian 
vector is an isotropically distributed unit vecfo/0. 

Assumption 1. The k*^ column of H, i.e., h^, has the 
same distribution as ax, where x is a normalized complex 
Gaussian vector, and a is a scaling factor 

Assumption [T] is valid for most practical MIMO systems. 
Consider the singular value decomposition (S VD) of Hk.k = 
UAV^, where U and V are unitary matrices containing 
the left and right singular vectors, respectively, and A is a 
diagonal matrix containing the singular values. It is known 
that the columns of U and V have the same distribution 
as a normalized complex Gaussian vector ITSl . In the case 
of single-user precoding, it is optimal to set to be 
proportional to vi, the right singular vector corresponding 
to the maximum singular value Ai lfT4l . As a result, is 
proportional to AiUi, where Ui is the left singular vector 
corresponding to the maximum singular value. Hence, As- 
sumption[T]is valid. It can also be shown that the assumption 
is valid when linear multiuser precoding is deployed. 

Assumption 2. The entries o/H(_i.) are independent com- 
plex Gaussian random variables. 

Assumption|2]is in general valid. Indeed, it is obvious that 

' A unit vector is said to be isotropically distributed if it is equally likely 
to point in any direction in the complex space. In other words, the vector 
is uniformly distributed on a complex unit sphere. 



the columns of 'ii-[-k) ^6 independent, as they are related 
to channel matrices on different links. Moreover, for a given 
tj, we have hj ~ CA/'(0,afej||tj||2l), thus implying that 
h.j has independent entries. Being the transmit precoding 
vector of the j*-^ primary link, tj is typically a function of 
Hj j and is independent of H^ j. This justifies the validity 
of Assumption 12] 

By Assumption|2] we can write 'ii-(-k) ^{-k) = HA, 
where 

A = diag(^Q;fea||ti||2, • ■ • , y/ak,k-i\\tk-i\\2, 

^Q;fc,fc+l||tfc+i||2, . . . , ^Jak,K\\tK\\■2), 

and H is an A^^ x {K— 1) matrix with independent standard 
complex Gaussian entries. Let H = UAV^ be the SVD of 
H. It is known that U and V are isotropically distributed 
unitary matricefl and that U, V, A are independent lfT3l . 

Proposition 1. has the same distribution as a normalized 
complex Gaussian vector as long as the random Hermitian 
matrix is a unitarily invariant matri^ 

Proof: Being a unitarily invariant matrix, can be 
decomposed as W/; = SDS^, where S is an isotropically 
distributed unitary matrix independent of the diagonal matrix 
D Us). Thus, Vk = ,,^^^Hi^'',,2 - Since the distribution of 

||SDS nk\\2 

is rotationally invariant and is independent of S, S^hfc 
has the same distribution as hfe. It follows that r^, has the 
same distribution as 

SDhfc g Dhfc 
||SDhfc||2 IIDhfclli' 

Upon conditioning on h^; and D, ^^^^''^^^ is a deterministic 
unit vector Since multiplying any deterministic unit vector 
by an isotropically distributed unitary matrix results in an 
isotropically distributed unit vector [15|, the unit vector 
is isotropically distributed for the given hfe and D. Since 
this holds for any realizations of and D, it follows that 
Yk is isotropically distributed, and therefore has the same 
distribution as a normalized complex Gaussian vector ■ 

Corollary 1. Vk has the same distribution as a normalized 
complex Gaussian Gaussian vector when a MF, ZF, or 
MMSE receiver is deployed. 

Proof: To prove Corollary [T] all we need is to show that 
Wfc is unitarily invariant for MF, ZF, and MMSE receivers. 
As this is trivial in the MF case, where — I, we will 
focus on the cases with ZF and MMSE receivers in the 
following. 

For ZF receivers, substituting H(-fc) — HA to (|6|, 
we_have Wj, I - UA[A^A)-i A^U^ = U(I - 
A(A^A)~^A^)U^. Since U is an isotropically distributed 

-A unitary matrix is said to be isotropically distributed if its probability 
density is unchanged when premultiphed by a deterministic unitary matrix. 

^A random Hermitian matrix W is called unitarily invariant if the joint 
distribution of its entries equals that of GWG^ for any unitary matrix 
G independent of W. 
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unitary matrix, it has the same distribution as GU for any 
unitary matrix G that is independent of U. Therefore, Wk 
is unitarily invariant. 

In the case of MMSE, ^ can be written as = 
U(T + Noiy^ij", where T = AV^A^VA is a Hermi- 
tian matrix. Since U is an isotropically distributed unitary 
matrix, Wfc is unitarily invariant. 



utiHzation of the spectrum, we can exploit the distribution 
of rk and consider the probabilistic interference constraints 

Pr I |rf Hfc,sts r <—}>l-6k Vfe = 1, . . . , i^. 

[' ak,s J 

(10) 

To tackle the constraints in (fTOl l. we need the following 
lemma: 



III. Optimal Beamforming as Homogeneous 
QCQP 

In this section, we show that the optimal SU beamforming 
problems that arise in the three scenarios discussed in 
SectionlUcan all be formulated as quadratically constrained 
quadratic programming (QCQP) problems. To begin, let us 
simplify Problem (|2|i by exploiting the properties of its 
optimal receive beamforming solution r^. Observe that the 
variable rs appears only in the objective function of Problem 
(|2]i. Thus, for any given ts, the optimal rs that maximizes 
75 is simply an MMSE receiver ||T6| given by 



r^(t5) = I3s^-'lis,sts, 



(8) 



where $ 
l|rj(ti 



||*-iHs,sts||2 

^ 1. 



s k ^ A'ol, and 
is a normalization factor that ensures 



Upon substituting r^(t5) into 75 (see ([T]i), we obtain 
75 - a5,5tfHf;s*-iHs,sts = tf Ats, where A = 
as,sHs s^^^^s.s- In particular, we can eliminate the 
variable rs from Problem Q and replace the objective 
function with the quadratic form tg Ats- 



A. Homogeneous QCQP Formulation in Scenario 1 

Now, recall that the secondary transmitter has perfect 
knowledge of H^^^r^ in Scenario 1. With the optimal 
given in Problem ^ becomes 



max tf Ats (9a) 
tfQfets<l Vfc = l,...,i^, (9b) 



ts 
s.t 



tfts < Ps, 



(9c) 



where = ^^^H^ grkrj^Hk,s are Hermitian positive 
semidefinite matrices. Problem (|9l) is a homogeneous QCQP, 
where both the objective function and inequality constraints 
are quadratic without linear terms. 

In subsequent subsections, we will show that similar ho- 
mogeneous QCQP problems can be formulated for Scenarios 
2 and 3, with in ( |9bl ) replaced by some suitable matrices 
and Q|, respectively. 

B. Homogeneous QCQP Formulation in Scenario 2 

In Scenario 2, the realization of rfc is not known to the 
secondary transmitter. In order to have a more efficient 



Lemma 1. Let r G C" be a normalized complex Gaussian 
vector, i.e., r = j^^^ where z ^ CJV{0,1) is a standard 
complex Gaussian vector. Let u e C" be an arbitrary vector, 
and let C > 0, (5 G (0, 1) be arbitrary scalars. Then, we have 



Pr{|r^up < C} > 1 



C 



1-6— 



Proof: Since the distribution of r is rotationally invari- 
ant, we may assume without loss of generality that u = 
[||u||2, 0, . . . , 0] . Then, we have 

Pr{|r^u|' <C} =Pr{||u||M[r]i|' <C} 



PrS (||u 



l2-C)|Ni|'< 



cEI 

1=2 



where denotes the i entry of a vector. Note that 
^i=2 '■^^ same distribution as ^X2{n-i)' where 

is the standard real chi-square random variable with d 
degrees of freedom. Moreover, it is independent of | [z] 1 1 , 
which has the same distribution as ^X2- Hence, we have 

Pr {|r^u|' < C} = Pr {(||u||^ - Qxl < Cx^(„-i)} 



pj 4n-l)/2(^-l) ^ ||U 



Now, let 

Ia{a,b) 



a+b-l 

E 

J=a. 



j!(a + 5-l-j)! ^ ^ 



(11) 

be the regularized incomplete beta function. It is known that 
the random variable 

X^(„-i)/2(n-l) 



F, 



2(n-l).2 



xi/2 



follows the so-called i^-distribution with (2(n — 1),2) de- 
grees of freedom, whose cumulative distribution function 
(CDF) is given by Pr{F2(„_i) 9 < x} = I 2(,,-i)x {n - 

-l)x + 2 



(n — l)x 



1,1). Upon working out the summation in (flTI) . we obtain 

Pr{^2(„-1),2 < x} = 



{n - l)a; + 1 
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Thus, we conclude that 

Pr||r^u|^ < C} > 1 - ^ 



< S 



|u||i/C-l 
l|u||i/C 



< S 



|U||2< 



c 



i-s— 

which completes the proof. ■ 
By Lemma [T] we can rewrite the chance constraints (fTOl l 

as 

1 



|Hfe,sts||^ < 



1-5 



k 



\fk = l,...,K, 



which of course are equivalent to 

tfQ^s<l Vfc = l, 

where 



(12) 



1 _ srl/{Wfc-l) 
^ "fc ^ xtH tt 



In particular, the optimal beamforming vector ts in Scenario 
2 can be found by solving the QCQP problem ( l9al ). (fT2]) and 
I©. 

As mentioned earlier, when 6k — 0, the chance constraints 
(ITOl l are equivalent to the worst-case constraints ((Sj. In this 
case, we have = ^lif^glik,s- 



C. Homogeneous QCQP Formulation and Closed-Form So- 
lution in Scenario 3 

1) Homogeneous QCQP: In Scenario 3, both and 
Hfc 5 are unknown to the secondary transmitter Similar to 
Scenario 2, we consider the following probabilistic interfer- 
ence constraints: 



Pr (|rfHfc,sts| 



< 



or equivalently. 



tfQlts<l Vfc = l,...,/^, 



(15) 



where Q| = log ^I. Thus, the optimal beamforming 
vector ts in Scenario 3 can be found by solving the QCQP 
problem dgali, ( fTsT i and ( l9c] i. 

We would like to emphasize that if the worst-case in- 
terference temperature constraints (O are to be satisfied in 
Scenario 3, i.e., 5k = 0, then the only feasible solution is 
ts = 0. That is, the secondary transmitter can never transmit 
under the overly stringent constraints. On the other hand, by 
allowing a small outage probability 5k, the secondary link 
can transmit on the same spectrum of the primary system. 

2) Closed-Form Solution: Here, we show that the op- 
timal ts in Scenario 3 can be found very efficiently by a 
simple eigenvalue-eigenvector computation. Indeed, observe 
that the constraints (fT4l i and i9c[ can be combined to yield 
the single constraint 



tsWl < A 



log ^ ' ■ ■ ■ ' aK,s log 



S.max 



Thus, the optimal beamforming problem in Scenario 3, 
which is given by 



max tgAtg 
ts 



s.t. 



|t5lli< A, 



(16) 



> 1-4 yk^l,...,K 
(13) 

Since r^^ and H^. g are independent, upon conditioning 
on Tfc, we see that r^H^ gtg is a complex Gaussian 
random variable with mean and variance Ijtslli- Note 
that the conditional distribution of r^VLk.s^s is inde- 
pendent of Tfe. It follows that unconditionally, we have 
r|^Hfe,sts ^ CA/'(0, lltslli)- In particular, the random 
variable jr^H^^stsp follows an exponential distribution 



is reduced to the problem of finding the largest eigenvalue of 
A and its associated eigenvector. Specifically, let v be the 
eigenvector of A corresponding to the largest eigenvalue. 
Then, the optimal solution to ( fTSI l is simply t*g = \/Av. 
Note that there is no need to solve any QCQP in this case. 



IV. SDP Relaxation 

We have shown in the last section that the optimal beam- 
forming solution can be efficiently obtained by a simple 
eigenvalue-eigenvector computation in Scenario 3. However, 
to obtain the optimal solutions for Scenarios 1 and 2, 
homogeneous QCQP problems of the following form have 
to be solved: 



max 
ts 



tfAts 



(17a) 



s.t. tfQfcts<l Vfc = l,...,i^, (17b) 



tf ts < P.. 



S.max • 



(17c) 



with parameter 



Its II 



Note that this result holds as long 



as Vk is a unit-length vector, regardless of its distribution. 
Since 



Pr r^H.sts < 



e-k 



1— cxp 



it follows that the chance constraints in ( [13] ) can be written 



as 



CSII2 



< 



ak,s log j- 



yk=l,...,K, (14) 



Here, is equal to and Q^, in Scenarios 1 and 
2, respectively. Unfortunately, since Problem ([TtI i involves 
maximizing a convex function over an intersection of K 
ellipsoids, it is NP-hard in general IfTTl . In this section, we 
show how Problem ( [TtI i can be tackled using semidefinite 
programming (SDP) relaxation methods. 

To begin, observe that t;5Qt5 = tr(QX) for any matrix 
Q, where X = tgt^ is a rank one Hermitian positive 
semidefinite matrix. Thus, by relaxing the rank constraint 
lank(X) = 1, we obtain the following SDP relaxation of 
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Problem ( fTTj i: 



max tr(AX) 

X)-0 



(18a) 



mm 

vi,---,yK ,yK+i 



solution to M7^ can be constructed from an optimal solution 
of (II8I 1 in polynomial time. 

Proposition|2]establishes the existence of polynomial-time 
algorithms for constructing optimal solutions to (fTTI i when 
K <2.1vi subsections IV-AI and IV-BI we describe two such 
algorithms — one for the case where K = 1, and the other 
for the case where K ~ 2. Both of them are based on the 
following decomposition theorem of Huang et al. lfT9l . Il20l : 

Theorem 1. Suppose that Z is a Hermitian positive semidef- 
inite matrix of rank R, and A and B are two given Hermitian 
matrices. Then, there is a rank-one decomposition of Z, 
fc=i namely, Z = ^i-z^, such that z^Az^ ~ '^^d 

yk>0 Vfc = 1, . . . , A' + l.(19c) Bz,. = ^^^^ for all r ^ 1,2, R. Moreover, such a 



s.t. tr(QfcX) < 1 \fk^l,...,K, (18b) 
tr(X) < Ps ,7nax ■ 

(18c) 

The dual of (fTSl l is given by 

K 



.max VK+i, (19a) 



fc=i 

K 



S.t. J2 yk^k + VK+il A to, (19b) 



It is known that SDP problems are convex and can be solved 
in polynomial time using standard interior-point methods 
lITSl . Moreover, if we can find a rank-one optimal solution 
to (fTSl l, then we can extract from it an optimal solution 
to the original QCQP problem ( fTTb . In this case, there 
is no gap between the optimal value of iT% and that of 
(fTsT l, and Problems iT% and ( fTsT i are equivalent. Of course, 
the SDP relaxation dTsl ) is in general not equivalent to 
the QCQP problem dTTI i. as we have discarded the rank 
constraint rank(X) = 1. In the next two sections, we 
will discuss how to recover a rank-one solution from an 
optimal solution to ( fTSl l. Before proceeding, however, let us 
introduce the following lemma, which will be useful for our 
later discussions. 

Lemma 2. Both ( flST l and its dual (fT9] l satisfy the Slater 
condition, i.e., they are strictly feasible. 

We omit the proof here due to page limit. 

Remark 1. Since both (fTsl l and (fT9] l are strictly fea- 
sible, a pair of primal and dual feasible solutions 
(X*; {ul,. . . ^Vk+i)) to ( fTsl l and ( fT9] l is optimal if and only 
if the following complementary conditions hold: 

tr (x.* ylQk + yx+il - A^ ^ = 0, (20a) 

yl (tr(QfcX*) - 1) = Vfc = 1, . . .,K,(20h) 
y*K+i (tr(X*) - Ps,,naT) - 0. (20c) 

V. Optimal Rank-One Solution when K <2 

As it turns out, when there are no more than two primary 
links (i.e., when K < 2), there is no gap between the 
optimal value of the SDP relaxation ( fTSl l and that of the 
original QCQP problem ( [TtI i. Moreover, a rank-one optimal 
solution to ( fTsl l. and hence an optimal solution to ( fTTj l. 
can be found in polynomial time. Specifically, we have 
the following proposition, which follows directly from the 
results of Huang et al. Qll, l20l : 

Proposition 2. The homogeneous QCQP problem (I17l l can 
be solved exactly in polynomial time when the number of 
primary links K is at most 2. In particular, an optimal 



decomposition can be found in polynomial time. 

We refer the interested readers to |f20l for the proof. To be 
self-contained, the algorithm for computing the decomposi- 
tion guaranteed by Theorem [T] which runs in polynomial 
time, is given in Algorithm [T] 



A. Optimal Rank-One Solution when K = 1 

When K = \, there are two quadratic constraints in ( fTTb . 
namely, t^Qit^ < 1 and t^ts < Ps,max. Let X* be an 
arbitrary optimal solution to ( fTSI ), which can be obtained in 
polynomial time by standard interior-point algorithms. We 
now show how to construct a rank-one solution to ( fTSi ) from 
X*. By Theorem [U we can find a rank-one decomposition 

X* = J2r=l ^rtr SUCh that 



and tft.^^ 
Vr = l,...,i?, 



(21) 



where R is the rank of X*. Now, choose an arbitrary t^, 
say ti, and let t = v^ti. From ( f2TT i and the fact that X* 
is a feasible for ( fTSi ), we have 



t^Qit = tr(Qitt^) = tr(QiX*) < 1 



and i"i = tr(tt'') = tr(X*) < Ps. 



(22) 



Thus, we conclude from ( f22] i that tt^ is a rank-one feasible 
solution to (fTsl l. 

Now, from the complementary conditions in (f20l i. we have 

tr (X* {ylQ,+y;i-A)) 
= ^tr(t,tf (2/tQi+2/*I- A)) =0. 



Since yJQi + t/^ - A >z by (fT9b] i, we have 
tr(M^(yi*Qi+y2*I-A)) = for r = 1,...,R. 
This implies that tr (it" (ylQi + y^I - A)) = 0. In 
a similar fashion, since tr(Qitt^) = tr(QiX*) and 
tr(tt^) = tr(X*), we have yl (tr(Qitt^) - l) = 
and 2/2 (tr(tt-^) - Ps.rnax) = 0. This, together with the 
feasibility conditions in (f22] i. leads to the conclusion that 
the rank-one matrix it" is an optimal solution to (fTSl l, and 
t is an optimal solution to ( fTTI i. 
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Algorithm 1 Algorithm for computing the decomposition 
guaranteed by Theorem [T] 

Input: Hermitian matrices A and B, and Hermitian positive 

semidefinite matrix Z with R = rank(Z). 
Output: Z " 



S^Li ZrZ^, a rank-one decomposition of 



23: 
24: 

25: 
26: 
27: 
28: 
29: 



_ tr(AZ) 



_ tr(BZ) 



R 



for 



rank-one 



,.= 1 PrPr 



decomposition 
using. 



Z such that z^Az^ 
r = l,...,i?. 
Compute an arbitrary 
Pi,P2, . . . , Pfl such that Z 
for example, Cholesky factorization. 
Let r = 1. 
repeat 

if (pf Ap. - H(|^)(pf Ap, - > for all 

j = r + 1, . . . , R then 

Qr := Pr- 

else 

Let Z e 7' + 1, . . . , i? be such that (p^Ap,. - 

H(^)(pPAp,-ii(^)<0. 

Determine 7 such that (Pr +7P()^-^(Pr +7Pi) — 

^(1 + 7^). 

:= ^7==, and set p; 

end if 

if r = i? 1 then 



Pi- 



le 

17: 
18: 

19: 
20: 
21: 

22: 



end if 

r := r + 1. 

until r = i? — 1 {Comment: We have now found a 
decomposition Z = 2^=1 li-q^ such that q^Aq^ = 

tr(AZ) 1 
R 

Let r = 1. 



repeat 

if (q?Bq.-H(|^)(qfBq, 



tr(BZ) 



) > for all 



3 = r + 
Zr : = 

else 

Let / 

tr(BZ 



, i? then 



Note that t can be computed in polynomial time, as 
both X* and the decomposition X* = J2r=i *i-^r^ 
computed in polynomial time. 

B. Optimal Rank-One Solution when K — 2 

When K = 2, there are three quadratic constraints in 
(fTTl i. namely, tfQitg < 1, tfQats < 1 and tfts < 
Ps.max- In the following, we show how to construct a rank- 
one solution from an arbitrary optimal solution X* to (fTsT l 
in two different cases. 

1) At least one constraint in (22) is non-binding 
at optimality: Without loss of generality, suppose that 
tr(QiX*) < 1 is the non-binding constraint, while 
tr(Q2X*) < 1 and tr(X*) < Ps,max can be either binding 
or non-binding. Due to the complementary conditions ( I20bl i. 
we must have yl = 0. Now, construct, via Theorem [T] 
a rank-one decomposition X* = J2r=i ^r^r ^^^^ that 
t^Q^t, = and t«t,, = for r = 1, . . . , i?. 

Since tr(QiX*) = X^^Li^^Qitr, there must exist an 
r e such that t^Qit,. < Without 

loss of generality, assume that r — 1 and let t = \/i?ti. 
By the same argument as in the preceding subsection, the 
following feasible conditions and complementary conditions 
hold: 

tr(Q2tt^) < 1, tr(tP) < Ps,,nax: 

tr (tF (2/*Qi + 2/2* Q2 + 2/3*1 - A)) = 0, 



yl (tr(Q2tP) - 1) = 0, yl (tr(tt«) - 



S.max , 



= 0. 



Moreover, we have tr(Qitt^) = R ■ tr(Qititf ) < 
tr(QiX*) < 1 and y* (tr(Qitt^) - l) = since y* = 0. 
Hence, tt^ is an optimal rank-one solution to (fTSl l and t is 
an optimal solution to ([TtT i. 

2) All constraints in t^f^ binding at optimality: 
When all constraints are binding, we have tr(QiX*) = 
tr(Q2X*) = 1 and tr(X*) = Ps max- Then, we have 
tr((Qi - Q2)X*) = tr((Q2 - ^^I)X*) = 0. In 
this case, we construct, again via Theorem [T] a rank-one 



e r + 1, . . . , £ be such that (qf Aq^ - decomposition X* = VfL, Kt" such that 



R 



)(qfBq,-™)<0, 



Compute the arguments qi := arg(q^Aq;) and 
a2 ■= arg(q^Bqi) and the modulus 70 = 
|q^Bq;|, and determine 7 such that (q^Bq^ — 

tr(BZ1\ 2 1 o • / \ I Ht3 tr(BZ) 

-4^)7^+270 sm(Q:2-Q!i)7+qfBq; ^-j^ = 

0. 

Set w = je 



tf (Qi-Q2)t, 



and 



,i(ai+7r/2) 



Q2 



P 



S,max 



It, 



tr((Qi -Q2)X* 
R 



tr Q 



(23) 



Ps., 



I X^ 



R 



end if 

if r = i? 

Zi? := q;- 
end if 
until r = i? 1 



Vl+T 

1 then 



, and set q; 



(24) 

for r ~ 1, . . . ,R. Since tr(QiX*) = 1, there must exist an 
r e {1, . . . , R}, say r = 1, such that tr(Qititi'^) = s > 0. 
Let t = and consequently we have tr(Qitt^) = 1. 
This, together with (|23]l and (EHi, leads to tr(Q2tt^) = 1 
and tr(tt^) = Ps.max- Hence, tt^ is a feasible solution, 
and the complementary conditions in ( |20bt are satisfied. 
Furthermore, since {ylQi + 2/2Q2 + Vs^ ^ A) ^ 0, we see 
that tr (tt^ (ylQi + J/2Q2 + 2/3! - A)) = 0. Hence, we 
conclude that tt^ is an optimal rank-one solution to ( fTsl l. 
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and t is an optimal solution to the QCQP problem ( flTl l. 
Moreover, t can be computed in polynomial time. 



can decompose 



as Qfc = Pkj(^^kj^H 



VI. Rank-One Solution when K > 3 

When K > 3, there may not exist any rank-one optimal 
solution to the SDP ([TSll. Moreover, the QCQP problem ( [TtI i 
is NP-hard in general, and hence it is unlikely that we can 
extract, in polynomial time, an optimal solution to it from 
an optimal solution to the SDP (fTsl l. However, as we shall 
see, we can generate a provably near-optimal solution to the 
QCQP problem ([TtI i from an optimal solution to the SDP 
(fTsT l using a very simple randomized procedure. 

To begin, let X* be an optimal solution to the SDP 

-I, so that constraint ( |9c] i is 



Ps., 



(fTST i. Define Q^+i 
equivalent to tg QK+its < 1. Consider the randomized 
procedure outlined in Algorithm |2] for generating a feasible 
solution to ([TtT i from X* . Algorithm |2] can be viewed as 
a generalization of the procedure developed by Nemirovski 
et al. iflTl for handling real homogeneous QCQP problems. 
Our goal now is to show that Algorithm |2] indeed returns a 
feasible solution to ST% . In fact, we will prove in Theorem |2] 
that the solution returned by Algorithm|2]is not only feasible 
to ( [TtI i. but is also likely to be a good one, in the sense that it 
has an objective value that is close to the optimal value of the 
QCQP problem il7\ . We note that such a phenomenon can 
also be observed from our simulations, as will be explained 
in the next section. 

Algorithm 2 Generate a feasible solution to il7\ from an 
optimal solution X* to ( fTSl l 
Input: An optimal solution X* to the SDP (fTsT l. 
Output: A feasible solution t to (fTTI i. 
1: Decompose X* = A^A, where A e cMsxMs^ 
A = AAA^ and Qk = AQ^A^, where k = 
1, ...,/£+ 1. It can be shown that tr(A) tr(AX*) 
and tr(Qfe) = tr(QfeX*) < 1. 
2: Find an unitary matrix U that diagonalizes A, i.e., A = 

IJ-^AU is a diagonal matrix. Set Qfc = U-^QfcU. 
3: Let <^ be an A/5 x 1 random vector whose entries 
are independently and uniformly distributed on the unit 
circle in the complex plane. In other words, we have 
= c-'^*' where 9i is uniformly distributed between 
and 2tt. 
4: Return t 



A^U^ as the solution. 



Before we introduce and prove Theorem IH let us note the 
following facts: 

Fact 1. I \19i There exists an optimal solution to Problem 
(II8I 1 with rank R < ^/K + 1, where K +1 is the number of 
quadratic constraints. Moreover, such an optimal solution 
can be found in polynomial time. 



Fact 2. We have rank(Qfc) 
min{Vir + 1, Ms} /or fc = 1, . . 



= rank(Qfc) 

.K + 1. In particular, we 



In order to study the quality of the solution returned by 
Algorithmic we need the following lemmata: 

Lemma 3. Let a > be given. Consider the events 

K fi 
/c=lj = l 

where f '^^ is obtained from the rank-one decomposition of 
Qfc (see Fact|2]). Then, we have 

Pr I max ^^Qfc^ > al < Vy{A\. 

[l<fe<A'+l J 

Proof: If A does not take place, then we have 
^Hffej(-fA;j)if^ < aWi^^W for all fc = 1,...,A' + 1 and 
i ~ 1, . . . , fi. This implies that 



yk = l,...,K + 1. (25) 

Note that -£^=1 W^'^'Wl = tr(Qfe) = tr(Qfc) < 1. Hence, 
dZSl l implies that maxi<:k<K+i C'^QfeC ct when A does 
not take place. This completes the proof. ■ 

Lemma 4. (Hoeffding's Inequality, Complex Version) Let 
Xi, . . . , Xn be independent complex-valued random vari- 
ables with KXi = and < a.i for i = 1, . . . ,n. Then, 
for any /3 > 0, we have 



Pr 



1=1 



> P} <4: exp 



13' 



4 a 



where ||a||2 denotes the ^2-norm of the vector a = 
[ai, . . . , a„]'^. 

Proof: Let X[^ and Xj be the real and imaginary parts 
of X„ respectively. Then, we have EXf = EX/ = 0, 

l^i^l < o-i and \X.l I < Ui for i = 1, . . . ,n, and 



j:x. 

i=i 

n 



>f3 



> 



V2 



Pr 



> 



V2 



The desired result then follows from an application of the 
real version of the Hoeffding inequality j21;]. ■ 

We are now ready to present Theorem |2] It extends 
Nemirovski et al's result in iflTl . which is concerned with 
real homogeneous QCQP problems, to the case of complex 
homogeneous QCQP problems. 

Theorem 2. The vector t returned by Algorithm \2\ is well 
defined and is a feasible solution to Problem illi . Moreover, 
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for any a > 0, we have 

Prjt^At > ^tr(AX*)| > 1 - 4(X + 1)^ cxp , 

(26) 

where fi = min{\/^ + 1, Ms}- In particular, with proba- 
bility at least 1 — 4:{K + cxp (— j), the objective value 
of the solution returned by Algorithm |2] is at least ^ times 
the optimal value of the QCQP problem dTTl ). 

Proof: We first prove that t is well defined, i.e., 
maxi</j<A'+i C'^Qfc^ > 0. To see this, note that C^Q/c^ = 
^^U^AQfcA^U^ = t^Qfet, where t ^ A^U^ Since 
Qif_i_i >- 0, it follows that maxi<fe</f-|_i t^Q^t, and hence 
niaxi<fc<x+i C^QfcC' must be strictly larger than zero. 
Now, observe that 



t^Qfet = 



maxfc C^QfcC 



e^U^AQfcA^U^ 



1 



-A'^Qk^ < 1 



maxfe ^^QfcC 

for fc = 1, . . . , A' + 1. It follows that t is a feasible solution 
to (0. 

Next, we compute 

1 



t«At = 



1 



maxfc^^Qfc^ 



maxfc £,"Qk£. 
tr(A) = 



1 



maxfc^^Qfc^ 



tr(AX* 



where the second equaUty is due to the fact that A is a 
diagonal matrix and = 1 for i — 1, . . . , Ms- Hence, 

to prove the bound in (|26] |, it suffices to show that 



< 4(a: + i)^cxp (- 



Pr<! max ^^Qfe^ > a 

l<k<K+l 



Now, by Lemma [3] we have Pr |maxfe ^^Q^^ > a| 



(27) 



MA} < Efc,, 
and 



PilAkj}- Moreover, since E{[^]4f'=J]J = 
A/5, we have 



= Pr 



[fn.HI[f'^h|for* = l, 

Pr{|e^f''^-'"|> V^||f^--'||2} 

> V^llf'^-'llsj <4exp(-| 



Pr{Aj} 

Ms 



by Lemma m This establishes ( |27] | and hence the bound in 

Finally, the last statement in the theorem follows from the 
observation that tr(AX*) is an upper bound on the optimal 
value of the QCQP problem ( [TT] ). as (fTST l is a relaxation of 
(fTTl i. This completes the proof of Theorem [2] ■ 

Before leaving this section, we emphasize that the optimal 
beamforming vector can always be found efficiently in 
Scenario 3 through a matrix eigenvalue-eigenvector com- 
putation, regardless of the number of primary links. In 
Scenarios 1 and 2, however, the optimal solution can be 
obtained in polynomial time only when K is no larger than 
two. Otherwise, we can only find an approximate solution 
in polynomial time via Algorithmic Fortunately, as we will 



show in the next section, the approximate solution is nearly 
optimal most of the time. 

VII. Numerical Simulations 

In this section, the performance of the proposed algo- 
rithms are evaluated through simulations. Throughout this 
section, we assume that all stations are equipped with 4 an- 
tennae. The wireless fading channel is Rayleigh distributed, 
and path loss exponent equals 4. The length of the secondary 
link is 10 meters, and Ps,max is chosen in such a way that 
the average SNR received by each antenna at the secondary 
receiver is lOdB, if there is no interference. We also assume 
that all primary users transmit at power Ps.max- Likewise, 
the transmit beamforming vector of primary user k are set 
to be the dominant right singular vector of 'H.k,k- Meanwhile, 
the primary receivers use MMSE beamforming vectors, as 
given in Subsection II-C. Unless otherwise stated, Sk is set 
to 1% for all k in Scenarios 2 and 3. Each point in the 
figures is an average of 50000 independent simulation runs. 

A. K = 2 

We first investigate a network with one secondary link and 
two primary links. The primary Unks are placed such that 
the distances between the secondary transmitter and the two 
primary receivers are 15 and 13 meters, respectively, while 
the distances between the secondary receiver and the two 
primary transmitters are 12.4 and 12.7 meters. As discussed 
in previous sections, the optimal beamforming solution t*g 
can be found in polynomial time in this case. 

In Fig. [U the optimal SINR = (tJ)^At^ (in dB 
scale) is plotted against e^, when varies from to 10 
dB for all k. It is not surprising to see that 7^ increases with 
the increase of the tolerable interference at the primary 
receivers. Meanwhile, the more channel information at the 
secondary transmitter, the higher the SINR at the secondary 
receiver, especially when ej. is low. Noticeably, the SINR 
gap between the three scenarios narrows when increases. 
This is because when the primary users can tolerate higher 
interference levels, the secondary user can spend less effort 
in eliminating interference to the primary users. Hence, the 
advantage of knowing 5 and becomes less obvious. 

Fig. |2] illustrates the tradeoff between the optimal SINR 
7^ of the secondary link and the outage probability 5k of 
the primary links. It is not surprising that in both Scenarios 
2 and 3, the secondary link can achieve a higher SINR when 
the primary links can tolerate a higher outage probability. 



B. /V =4 

We now simulate a network with one secondary link and 
four primary links. The distance between the secondary 
transmitter and the four primary receivers are 20, 18, 15 
and 13 meters, while that between the secondary receiver 
and the four primary transmitters are 16, 14, 12.4 and 
13.2 meters, respectively. With four primary links, only 
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e/N„ (dB) 

Fig. 1. SINR at the secondary receiver vs. tf^/No when K = 2. 




Scenario 3 



<3 — Achieved SINR 
— Upper bound on the optimai SiNR 



e/N„ (dB) 

Fig. 3. SINR at the secondary receiver vs. efc/A^o when K = A. 



1£ 




— * — Scenario 2 
-< — Scenario 3 




Fig. 2. Tradeoff between SINR and outage probability Sf^ when t^/No 
bdB and K = 2. 



Fig. 4. Tradeoff between SINR and outage probability 5^. when e^/No 
5dB and K = 4. 



approximate solutions can be obtained in polynomial time 
in Scenarios 1 and 2, as discussed in Section |Vl] 

In Fig. |3] the randomized algorithm in Algorithm |2] is 
carried out to obtain the beamforming vector in Scenarios 
1 and 2. The optimal beamforming vector in Scenario 3 
is obtained through an eigenvalue-eigenvector computation. 
For comparison, we also plot the optimal value of the SDP 
relaxation dTsl l. which is an upper bound on the maximum 
achievable SINR. As the figure shows, the randomized 
algorithm performs very close to the optimum. The achieved 
SINR almost overlaps with the upper bound of the optimal 
SINR. Meanwhile, similar conclusions drawn from Fig. [T] 
also apply here. 

The tradeoff between 7^ and 5k in the four primary link 
case is illustrated in Fig. |4] Recall that in Scenario 3, the 
only feasible solution when Sk = is = 0. Fortunately, 
the achievable SINR in Scenario 3 increases rapidly with 6k 
as long as Sk > 0, as shown in both Fig. [T]and[3] 



C. A Grid Network with 9 Primary Links 

In this subsection, we consider a network with 9 primary 
links arranged in a 70-by-40 meter grid as shown in Fig. |5] 
The lengths of all link are equal to 10 meters. The secondary 
link is randomly placed in the area. In Fig. |6l SINR at the 
secondary receiver is plotted against Each point in the 
curve is an average of 20000 independent secondary-link 
placements. 

Again, the figure shows that Algorithm |2] performs very 
close to the optimum. The achieved SINR almost overlaps 
with the upper bound of its optimal value. With full CSI, 
Scenario 1 can achieve a much higher SINR than Scenarios 
2 and 3, especially when efc is small. The better performance, 
however, comes at a price. In practical systems where full 
CSI is not available, one has to resort to the schemes 
developed for Scenarios 2 and 3 to achieve the maximum 
SINR. 
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Fig. 5. Placement of 9 primary links. 
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-e — Achieved SINR 

—3^ — Upper bound on the optimal SINR 



e/N„ (dB) 

Fig. 6. SINR at the secondary receiver vs. e^/A^o when K = 9. 

VIII. Conclusions and Discussions 

In this paper, we considered optimal secondary-link 
beamforming in MIMO CR networks when the secondary 
transmitter has complete, partial, or no channel knowledge 
on the links to primary receivers. We proposed a unified 
homogeneous QCQP formulation for all three scenarios with 
either deterministic or probabilistic interference-temperature 
constraints. In Scenario 3, the QCQP problem reduces to a 
matrix eigenvalue-eigenvector computation problem, which 
can be solved very efficiently. For Scenarios 1 and 2, we 
approached the QCQP problem by SDP relaxation. Notably, 
the SDP relaxation admits no gap with the true optimal 
value when there are no more than two primary links. In this 
case, the optimal beamforming solution can be computed in 
polynomial time. When the number of primary users exceeds 
two, we proposed a randomized polynomial-time algorithm 
that can construct a provably near-optimal solution to the 
QCQP problem from an optimal solution to the SDP. 

The reader may notice that there is a gap between the 
theoretical performance of the randomized polynomial-time 



algorithm (Algorithm |2]i as established in Section |VI] and 
its practical performance as demonstrated in Section IVIII 
This can be attributed to the fact that the main theoretical 
result in Section [Vll namely Theorem |2] only provides 
a worst-case guarantee on the performance of Algorithm 
ID In other words, the guarantee is valid regardless of 
the distribution of the input data. However, in the setting 
of MIMO CR networks, the input data follow a specific 
probability distribution, and the worst-case instance may not 
arise too frequently. It would be interesting to see whether 
one can obtain better theoretical guarantees by performing 
a probabilistic analysis of the performance of Algorithm |2] 
(see 1221 and the references therein for related work). 

So far, we have considered one secondary link only. 
However, the proposed schemes can be easily extended to 
a multiple-secondary-link system with the aid of medium- 
access-control (MAC). Suppose that there is a narrowband 
busy-tone channel in addition to the data-transmission chan- 
nel. When a secondary link wishes to transmit a packet, 
it first senses the channel to see whether there is another 
secondary link transmitting. If not, it sends a short busy 
tone on the busy-tone channel to reserve the airtime. Other 
secondary links, upon hearing the busy tone, will keep silent 
during the airtime reserved by the transmitting link. Having 
successfully reserved the airtime, the link will then start to 
transmit its data packet on the data-transmission channel. 
In case more than one secondary transmitter sends busy 
tones at the same time, a collision has occurred on the 
busy-tone channel and the secondary transmitters will each 
wait for a random time period before attempting again. By 
doing so, it is guaranteed that there is only one secondary 
link transmitting data packets at a time, and the proposed 
optimal beamforming methods can be applied. For practical 
implementation, we can adopt the random-access protocols 
in IEEE 802.1 1 wireless local area networks (WLANs), such 
as RTS/CTS DCF, to coordinate the contention on the busy- 
tone channel. 

Note that multiple secondary links can also coexist with- 
out the aid of a MAC protocol by properly configuring their 
beamforming vectors, preferably in a distributed manner 
In this case, secondary links interfere with each other, and 
thus the optimal beamforming problem becomes much more 
challenging. We will address this problem in our future 
research. 
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